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AN  APPLICATION  OF  TENSOR  THEORY  TO  3-D  SHAPE  ANALYSIS^ 

NIRWAN  ANSARI 

Purdue  University 
School  of  Electrical 
West  Lafayette,  IN 
(317)  494-1740 

Abstract 

In  this  paper,  -we  present  two  new  approaches  to  address  the  shape  recognition  problem.  The 
first  concept  is  an  extension  of  2-0  landmark* based  shape  analysis  to  the  3-D  case.  Secondly,  we 
extend  . the  work  of  Cyganski  and  Orr  to  rigid  non-planar  objects.  Both  approaches  will  be 
mathematically  formulated,  however  no  experimental  results  will  be  shown  in  this  paper. 

I.  3-D  Landmark- baaed  Shape  Analysis 

Objects  can  be  represented  by  landmarks  or  tokens  In  3-D  space.  Landmarks  or  tokens  are 

points,  usually  defined  along  the  object  surface,  that  hare  important  shape  attributes.  These  attri¬ 
butes  are  chosen  a  priori  and  are  problem  dependent.  Example  of  landmarks  are  corners  and  points 
of  high  curvature.  We  shall  assume  throughout  this  paper  that  the  landmark  locations  have  been 
detected.  We  shall  address  the  problem  of  processing  the  landmarks  to  identify  the  object.  To  iden¬ 
tify  an  observed  /unknown)  object  from  a  library  of  objects,  we  find  the  necessary  transformation 
that  maps  landmarks  of  each  library  object  to  the  corresponding  landmarks  of  the  observed  object. 
It  is  difficult,  even  if  it  exists,  to  determine  this  transformation.  If  the  transformation  can  be  found, 
then  shape  features  based  on  the  transformation  can  be  extracted  and  used  for  object  identification. 
W«  shall  address  this  problem  similar  to  the  way  we  addressed  the  problem  of  2-D  shape  analysis 
(AnD85|.  For  2-D  shape  analysis  the  object  was  partitioned  into  triangles  and  the  mapping  used  was 
triangular  deformation  obtained  through  an  affine  transformation.  The  shape  features  in  the  2-D 
case  were  based  on  the  dilatations  obtained  from  the  transformation.  In  3-D  we  partition  the  object 
into  tetrahedra  and  use  the  tetrahedral  deformation,  i.e.  a  transformation  from  one  tetrahedron  to 
another. 
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1.1.  Quaaiconformal  mapping 

Tetrahedral  deformation  can  be  determined  by  a  3-D  affine  transformation.  The  3-D  affine 
transformation  is  a  class  of  quaaiconformal  mappings  in  3-D  which  are  defined  [Vaifil]  as  follows: 

A  difieomorphism,  Mi—* fl,  (H,  H  C  R3)  is  called  K-quasiconformal  if 

D(x)  <  K  <  oo  for  some  constant  K, 

D(x)  -  max(D(x)Jl(x))  , 


where 


-j-5j — t - 

UW! 


DW¬ 


DM 


lilslL 


min 

:*!«] 


;r(x)ei)3 


max 


J(x)  and  P  are  the  Jacobian  and  the  derivative  of  f  respectively.  Note  that  ,  . 


!f(x)ej  and 


min 


,  ,  .Jf(x)«!  correspond  to  the  largest  and  the  smallest  eigenvalue  of  the  Jacobian  matrix  of  f. 
Effx),I2(x),  and  D(x)  are  known  as  the  outer,  inner  and  maximal  dilatations  of  f,  and  are  related  to 
the  eigenvalues  of  the  Jacobian  matrix  of  f  by  the  above  equations.  Infinitesimally  f  maps  a  ball  to 
an  ellipsoid.  These  dilatations  at  each  x«i"l  are  various  combinations  of  the  rstios  of  the  semi-axes  of 
the  3-D  ellipsoid.  In  general,  it  is  difficult  (if  it  exists)  to  find  a  diffeomorphic  map  of  a  Library  object 
to  an  observed  object.  The  affine  transformation  from  one  tetrahedron  to  another  is  however 
uniquely  determined.  Since  the  3-D  affine  transformation  is  a  class  of  quaaiconformal  mappings,  we 
can  use  the  above  definitions  to  obtain  the  dilatations.  The  maximal  dilatation  derived  from  the 
affine  transformation  can  be  used  as  a  distortion  measure  between  the  two  tetrahedra  (they  are  simi¬ 
lar  if  the  maximal  dilatation  is  L).  Using  an  approach  similar  to  that  of  partitioning  a  2-D  object 
into  triangles,  we  partition  a  3-D  object  into  tetrahedra.  We  then  evaluate  the  maximal  dilatation 
and  the  scale  change  (magnitude  of  the  Jacobian)  of  each  tetrahedral  deformation  between  a  library 
object  and  the  observed  object,  [f  the  two  objects  are  similar,  for  each  tetrahedral  deformation 
between  the  two  objects,  the  maximal  dilatation  is  one  and  the  scaie  change  remains  constant. 
Hence,  the  standard  deviations  of  the  maximal  dilatation  and  the  scale  change  associated  with  all 
the  tetrahedral  deformations  between  the  two  objects  can  be  used  to  determine  how  similar  the  two 
objects  are.  The  smaller  the  deviations  are,  the  more  similar  the  two  objects  are. 

1.2.  Strain  tensors 

The  above  concept  can  similarly  be  formulated  by  tensor  analysis.  We  assume  that  the 
observed  object  is  deformed  from  a  library  ooject.  Denote  . .  (with  respect  to  reference  frame  X. 
x\  i— 1.2.3)  as  the  region  occupied  by  the  library  object,  and  fl  I  with  respect  to  reference  frame  X, 

L  Thi»  research  waa  partially  supported  by  the  Air  Force  Office  ot'  Scientific  Retearch  uoder  Grant  F4902O- 
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?,  i=  1,2,3)  as  that  of  the  observed  object.  We  assume  that  the  deformation  from  ft  to  fl  is 
homeomcrphic.  Assuming  both  X  and  X  are  orthogonal  cartesian,  the  square  of  the  arc  length  of  the 
library  and  observed  object  using  the  general  summation  rule,  are  respectively 
ds*  »  dx'dx1  and  ds'  =  dx'dx1  i=  1,2,3.  ^ 

Thus, ds‘ -  ds: »  2^|idxidxi  ij*l,2,3.  and  n-,-t  »  - ^j)- 

2  d xi  dx1 

The  tensor  which  la  symmetric,  is  called  the  strain  tensor  [Sok51j.  The  eigenvalues  of  are 
called  the  principal  strains,  and  the  corresponding  eigenvectors  are  orthogonal.  The  directions 
corresponding  to  the  principal  strains  are  called  the  principal  directions  (axes)  of  the  strain  tensor. 
At  any  given  point,  it  is  a  mapping  from  a  ball  to  an  ellipsoid.  The  ratio  of  each  semi-axis  of  the 
ellipsoid  to  the  diameter  of  the  ball  equals  Vl+2‘\i  i=  1,2,3,  where  X;’s  are  the  eigenvalues  of  fiK. 
Again,  the  deformation  from  a  library  object  to  the  observed  object  may  not  be  homeomorphic.  It  is 
difficult  to  determine  the  deformation  function  even  if  it  is  homeomorphic.  The  tetrahedral  deforma¬ 
tion  used  for  the  above  landmark- based  problem  can  however  be  uniquely  determined.  The  various 
dilatations  (inner,  outer,  maximal)  of  the  tetrahedral  deformation  can  also  be  easily  computed  by 
means  of  the  eigenvalues  of  the  strain  tensor  associated  with  the  tetrahedral  deformation.  The 
strain  tensor  analysis  which  provides  an  intuitive  analysis  of  deformation  thus  can  also  be  used  to 
solve  the  landmark-based  problem. 

IT.  Surface  tensors 

We  shall  develop  another  concept  for  object  identification  dealing  with  objects  that  are  charac¬ 
terised  by  surface  patches  rather  than  sets  of  landmarks.  Object  identification  becomes  a  problem  of 
identifying  surface  patches.  Let  0  with  coordinate  reference  frame  X  be  a  surface  occupied  by  a 
library  object,  and  fl  with  coordinate  reference  frame  X  be  a  surface  of  an  observed  object.  This 
notation  is  consistent  with  the  strain  tensor  notation  presented  above,  the  difference  being  that  ft  is 
a  surface  instead  of  a  set  of  landmarks.  We  further  assume  that  fl  and  fl  are  related  by  an  affine 
transformation,  i.e.  5?  -  arjx1— c  )  ?,  where  c‘  and  c  are  the  centers  of  gravity  of  Cl  and  f7.  Here,  we 

are  interested  in  finding  the  a/,  the  coefficients  of  the  affine  transformation.  We  can  thus  transform 
the  library  object  (surface)  with  respect  to  the  observed  object  (surface)  by  the  affine  transformation. 
The  similarity  between  the  library  object  (surface)  and  the  observed  object  (surface)  is  determined  by 
how  the  transformed  library  object  (surface)  is  correlated  with  the  observed  object  (surface).  We 
describe  how  to  determine  the  affine  coefficients  below.  We  form  a  surface  (moment)  tenaor  on  C1  and 
H  as  follows: 

S'1"  *  J^(x*-c‘)(xi-ci)  •  •  •  dx'dsrdx3  -  f  (x‘-c‘)(x‘-c>)  ■  •  ■  ln(x‘ rxr^)dxldrdx3 

/n(?-?)(5?-?)  •  •  •  di'dird?  -  /(?-?)(?-?)  •  •  •  yi'x^di'dirdx3 

/  ai(xi-cl)a(»(x,-ci)  •  •  •  |  J’dx'dbrdr3  =*  a^a)  *  •  • 
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Thus,  S'1"-  is  a  relative  tqnsor  of  weight  -l.  Following  the  above  steps  we  form  various  orders  of  sur¬ 
face  tensors  for  both  Cl  and  fl,  we  then  contract  them  to  unit  rank  absolute  tensors  similar  to 
(CyOSSl.  Three  pairs  of  such  unit  rank  absolute  tensors  of  ft  and  ft.  say,  u‘,  IT,  v‘,  v,  w*,  w1,  are 
required  to  form  a  system  of  linear  equations: 

T  =»  a;uJ,  ~  *  af'v>,  w-  =  arw* 

such  that  we  can  obtain  ai. 
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